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BLACK-SCHOLES IN A CEV RANDOM ENVIRONMENT: A NEW APPROACH TO 

SMILE MODELLING 

ANTOINE JACQUIER AND PATRICK ROOME 


Abstract. Classical (Ito diffusions) stochastic volatility models are not able to capture the steepness of small- 
maturity implied volatility smiles. Jumps, in particular exponential Levy and affine models, which exhibit 
small-maturity exploding smiles, have historically been proposed to remedy this (see m for an overview). A 
recent breakthrough was made by Gatheral, Jaisson and Rosenbaum m, who proposed to replace the Brownian 
driver of the instantaneous volatility by a short-memory fractional Brownian motion, which is able to capture 
the short-maturity steepness while preserving path continuity. We suggest here a different route, randomising 
the Black-Scholes variance by a CEV-generated distribution, which allows us to modulate the rate of explosion 
(through the CEV exponent) of the implied volatility for small maturities. The range of rates includes behaviours 
similar to exponential Levy models and fractional stochastic volatility models. As a by-product, we make a 
conjecture on the small-maturity forward smile asymptotics of stochastic volatility models, in exact agreement 
with the results in m for the Heston model. 


1. Introduction 

We propose a simple model with continuous paths for stock prices that allows for small-maturity explosion 
of the implied volatility smile. It is indeed a well-documented fact on Equity markets (see for instance [29l 
Chapter 5]) that standard (ltd) stochastic models with continuous paths are not able to capture the observed 
steepness of the left wing of the smile when the maturity becomes small. To remedy this, several authors have 
suggested the addition of jumps, either in the form of an independent Levy process or within the more general 
framework of affine diffusions. Jumps (in the stock price dynamics) imply an explosive behaviour for the small- 
maturity smile and are better able to capture the observed steepness of the small-maturity implied volatility 
smile. In particular, Tankov |58j showed that, for exponential Levy models with Levy measure supported on 
the whole real line, the squared implied volatility smile explodes as o-^(fc) ^ —fc^/(2rlogr), as the maturity t 
tends to zero, where k represents the log-moneyness. 

Gatheral, Jaisson and Rosenbaum [30] have recently been revisiting stochastic volatility models, where the 
instantaneous variance process is driven by a fractional Brownian motion. They suggest that the Hurst exponent 
should not be used as an indicator of the historical memory of the volatility, but rather as an additional 
parameter to be calibrated to the volatility surface. Their study reveals that H G (0,1/2) (in fact H Ki 0.1 
in their calibration results), indicating short memory of the volatility, thereby contradicting decades of time 
series analyses. By considering a specific fractional uncorrelated volatility model, directly inspired by the 
fractional version of the Heston model [la |34], Guennoun, Jacquier and Roome |33) provide a theoretical 
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justification of this result. They show in particular that, when H G (0,1/2), the implied volatility explodes as 
cr‘^{k) ~ /T{H + 3/2) as t tends to zero (where ?/o is the initial instantaneous variance). 

In this paper we propose an alternative framework: we suppose that the stock price follows a standard Black- 
Scholes model; however the instantaneous variance, instead of being constant, is sampled from a continuous 
distribution. We first derive some general properties, interesting from a financial modelling point of view, and 
devote a particular attention to a particular case of it, where the variance is generated from independent CEV 
dynamics: Assume that interest rates and dividends are null, and let S denote the stock price process starting at 
S'o = 1, the solution to the stochastic differential equation dSr = Sr^/VdWr, for r > 0, where W is a standard 
Brownian motion. Here, V is a random variable, which we assume to be distributed as V ~ Yj, for some t > 0, 
where Y is the unique strong solution of the CEV dynamics dY^ = Yb > 0 where p € R, ^ > 0 and 

B is an independent Brownian motion (see Section 12.11 for precise statements). The main result of this paper 
(Theorem 12.31) is that the implied volatility generated from this model exhibits the following behaviour as the 
maturity r tends to zero: 


( 1 . 1 ) 


rlik) 




2(1-p) / fcV(l-p)t 

3 - 2p V 2r 
k'^^H 

r(logT)2’ 


l/(3-2p) 


ifp < 1, 
ifp= 1, 


k2 

. 2(2p-l)r|logr|’ ifp>l, 

for all fc ^ 0. Sampling the initial variance from the CEV process at time t induces different time scales 
for small-maturity spot smiles, thereby providing flexibility to match steep small-maturity smiles. For p > 1, 
the explosion rate is the same as exponential Levy models, and the case p < 1/2 mimics the explosion rate 
of fractional stochastic volatility models. The CEV exponent p therefore allows the user to modulate the 
short-maturity steepness of the smile. 

We are not claiming here that this model should come as a replacement of fractional stochastic volatility 
models or exponential Levy models, notably because its dynamic structure looks too simple at first sight. 
However, we believe it can act as an efficient building block for more involved models, in particular for stochastic 
volatility models with initial random distribution for the instantaneous variance. While we leave these extensions 
for future research, we shall highlight how our model comes naturally into play when pricing forward-start 
options in stochastic volatility models. In [41] the authors proved that the small-maturity forward implied 
volatility smile explodes in the Heston model when the remaining maturity (after the forward-start date) 
becomes small. This explosion rate corresponds precisely to the case p = 1/2 in (11.11) . This in particular 
shows that the key quantity determining the explosion rate is the (right tail of the) variance distribution at the 
forward-start date (here corresponding to t). 

The paper is structured as follows: in Sections 12.11 and 12.21 we introduce our model and relate it to other 
existing approaches. In Section [2.3l we use the moment generating function to derive extreme strike asymptotics 
(for some special cases) and show why this approach is not readily applicable for small and large-maturity 
asymptotics. Sections 12.41 and 12.51 detail the main results, namely the small and large-maturity asymptotics of 
option prices and the corresponding implied volatility. Section 12.61 provides numerics and Section 12.71 describes 
the relationship between our model and the pricing of forward-start options in stochastic volatility models. 
Section O also includes a conjecture on the small-maturity forward smile in stochastic volatility models. 
Finally, the proofs of the main results are gathered in Section |3 








BLACK-SCHOLES IN A CEV RANDOM ENVIRONMENT: A NEW APPROACH TO SMILE MODELLING 


3 


Notations: Throughout the paper, the ~ symbol means asymptotic equivalence, namely, the ratio of the 
left-hand side to the right-hand side tends to one. 


2. Model and main results 

2.1. Model description. We consider a filtered probability space (fl, (J^s)s>o, P) supporting a standard 
Brownian motion, and let (Zs)s>o denote the solution to the following stochastic differential equation: 


( 2 . 1 ) 


dZ, = - i Vds + VVdWs, Zo = 0, 


where V is some random variable, independent of the Brownian motion W. The process {Zs)s>o, in finance, 
corresponds to the logarithm of the underlying stock price. This is of course a simple example of a stochastic 
differential equation with random coefficients, existence and uniqueness of which were studied by Kohatsu-Higa, 
Leon and Nualart [43] and Albs, Leon and Nualart |2]. In the case where V is a discrete random variable, this 
model reduces to the mixture of distributions, analysed, in the Gaussian case by Brigo and Mercurio [TOl ITT] . 
In Appendix | b 1 we discuss how to hedge in this model. In a stochastic volatility model where the instantaneous 
variance process (Vt)t>o is uncorrelated with the asset price process, the mixing result by Romano and Touzi [54] 
implies that the price of a European option with maturity r is the same as the one evaluated from the SDE (12.11) 
with V = T~^ fj 14ds. As T tends to zero, the distribution of V approaches a Dirac Delta centred at the initial 
variance Vq . Asymptotics of the implied volatility are well known and weaknesses of classical stochastic volatility 
models are well documented [53]. Although such models fit into the framework of ( 12 . 11 ) . we will not consider 
them further in this paper. Define pathwise the process M by Ms := —^s -I- Wg and let {Ts)s>o be given by 
Ts ■= sV. Then T is an independent increasing time-change process and Z = M 7 -. In this way our model can 
be thought of as a random time change. Let now iV be a Levy process such that (e^'’)s>o is a (J^s)s>o-adapted 
martingale; define V := t~^ fj 14 ds where E is a positive and independent process, then )s>o is a classical 
time-changed exponential Levy process, and pricing vanilla options is standard [M Section 15.5]. However, as 
the maturity r tends to zero, V converges in distribution to a Dirac Delta, in which case asymptotics are well 
known [58] . 

The model (12.11) is also related to the Uncertain Volatility Model of Avellaneda and Paras [5] (see also [TU[351 
148] 1. in which the Black-Scholes volatility is allowed to evolve randomly within two bounds. In this framework, 
sub-and super-hedging strategies (corresponding to best and worst case scenarios) are usually derived via the 
Black-Scholes-Barenblatt equation, and Fouque and Ren m recently provided approximation results when the 
two bounds become close to each other. One can also, at least formally, look at ( 12 . 11 ) from the perspective of 
fractional stochastic volatility models, first proposed by Comte et al. in [14], and later developed and revived 
in [Ts] [28l |30l [33]. In these models, standard stochastic volatility models are generalised by replacing the 
Brownian motion driving the instantaneous volatility by a fractional Brownian motion. This preserves the 
martingale property of the stock price process, and allows, in the case of short memory (Hurst parameter H 
between 0 and 1/2) for short-maturity steep skew of the implied volatility smile. However, the Mandelbrot-van 
Ness representation [33] of the fractional Brownian motion reads 

'■* diu, 


:= 


1 


1 


dWs 


Jo (t-sp J-oo\it-sp i-s)\ 
for all t > 0, where 7 := 1/2 —H. This representation in particular indicates that, at time zero, the instantaneous 
variance, being driven by a fractional Brownian motion, incorporates some randomness (through the second 
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integral). Finally, we agree that, at first sight, randomising the variance may sound unconventional. As 
mentioned in the introduction, we see this model as a building block for more involved models, in particular 
stochastic volatility with random initial variance, the full study of which is the purpose of ongoing research. 
After all, market data only provides us with an initial value of the stock price, and the initial level of the 
variance is unknown, usually let as a parameter to calibrate. In this sense, it becomes fairly natural to leave 
the latter random. 

The framework constituted by the stochastic differential equation (EH) is a simple case of a diffusion in random 
environment. We refer the interested reader to the seminal paper by Papanicolaou and Varadhan |53j . the 
monographs by Komorowski et al. |44) . by Sznitman [57], and the lectures notes by Bolthausen and Sznitman [ 8 ] 
and by Zeitouni EOj. The classical set-up (say in is that of a given probability space (r2,A, Q) describing 
the random environment and a group of transformations (Ta:)a;gR‘ii jointly measurable in x € and w G 17 
(the transformation essentially indicates a translation of the environment in the x-direction). Consider now two 
functions 6 , a : 17 —>■ and define 

b{x,uj) := b{Tx{uj)) and a{x,uj) := a{Tx{uj)), for all x G £ 17. 

For each x £ and a; £ 17, we let Qx,lj denote the unique solution to the martingale problem starting at x 
and associated to the differential operator £“ := ^a{-,uj)A + 6 (-,w)V. The probability law is called the 
quenched law, and one can dehne the solution to the corresponding stochastic differential equation Qa;_(^-almost 
surely. The annealed law is the semi-product Qa, := Q x Qx,u:, and corresponds to averaging over the random 
environment. Most of the results in the literature, using the method of the environment viewed from the particle, 
however, impose Lipschitz continuity on the drift b{) and the diffusion coefficient a(), and uniform ellipticity 
of a(). These conditions clearly do not hold for (12.1|) . where b{-,uj) = —^uj and a(-,uj) = w, since w takes values 
in [0, c»). We shall leave more precise details and applications of random environment to future research. As far 
as we are aware, this framework has not been applied yet in the option pricing literature, the closest being the 
recent publication by Spiliopoulos [56] , who proves quenched (almost sure with respect to the environment) large 
deviations for a multi-scale diffusion (in a certain regime), assuming stationarity and ergodicity of the random 
environment. Interestingly, diffusions (specifically Ornstein-Uhlenbeck processes) in random environments have 
been justified by Follmer and Schweizer [23] and by Horst [35] (see also [52]) as the continuous-time limit of 
discrete-time dynamics for the asset price defined as the outcome of equilibria determined by the interaction of 
informed and noise traders. 

2.1.1. Moment generating function. In [211 EH [37], the authors used the theory of large deviations, and in 
particular the Gartner-Ellis theorem, to prove small-and large-maturity behaviours of the implied volatility in 
the Heston model and more generally (in m) for affine stochastic volatility models. This approach relies solely 
on the knowledge of the cumulant generating function of the underlying stock price, and its rescaled limiting 
behaviour. For any r > 0, let A^{u,t) logE(e“'^^) denote the cumulant generating function of Z^, dehned 

on the effective domain := {m £ R. : |A'^(it,r)| < oo}; similarly denote A^(m) = logE(e“^), whenever it is 
well defined. A direct application of the tower property for expectations yields 

(2.2) A^{u,t) = A^ ^ 'a('a ^ 1)'^ ^ ^ ^ ^ 

Unfortunately, the cumulant generating function of V is not available in closed-form in general. In Section 12.31 
below, we shall see some examples where such a closed-form solution is available, and where direct computations 
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are therefore possible. We note in passing that this simple representation allows, at least in principle, for 
straightforward (numerical) computations of the slopes of the wings of the implied volatility using Roger Lee’s 
Moment Formula [3S] (see also Section r2.3.2l) . The latter are indeed given directly by the boundaries (in R) of 
the effective domain of A^. Note further that the model could be seen as a time-changed Brownian motion 
(with drift); the representation (12.21) clearly rules out the case where Z is a simple exponential Levy process (in 
which case would be linear in r). In view of Roger Lee’s formula, this also implies that, contrary to 

the Levy case, the slopes of the implied volatility wings are not constant over time in our model. 


2.2. CEV randomisation. As mentioned above, this paper is a first step towards the introduction of ‘random 
environment’ into the realm of option pricing, and we believe that, seeing it ‘at work’ through a specific, yet 
non-trivial, example, will speak for its potential prowess. We assume from now on that V corresponds to the 
distribution of the random variable generated, at some time t, by the solution to the CEV stochastic differential 
equation dYu = ^Y^dBu^ Lb = yo > 0 where p € R, > 0 and R is a standard Brownian motion, independent 
of W. The CEV process [51112] is the unique strong solution to this stochastic differential equation, up to the 
stopping time := infK>o{Wi = 0}. The behaviour of the process after depends on the value of p, and shall 
be discussed below. We let T{n;x) := r{n)~^ denote the normalised lower incomplete Gamma 

function, and m^ := ^{Yt = 0) = P(V = 0) represent the mass at the origin. Straightforward computations show 
that, whenever the origin is an absorbing boundary, the density Cp(y) = IP(^i S dy)/dy is norm decreasing and 
{p defined in (12.41) 1 


(2.3) 


mt = 1 - r 



vT-^^ \ 

2e{i-p)Hj 


>0; 


otherwise mt = 0 and the density Cp is norm preserving. When p G [1/2,1), the origin is naturally absorbing. 
When p > 1, the process Y never hits zero P-almost surely. Finally, when p < 1/2, the origin is an attainable 
boundary, and can be chosen to be either absorbing or reflecting. Absorption is compulsory if Y is required to 
be a martingale [36l Chapter III, Lemma 3.6]. Here it is only used as a building block for the instantaneous 
variance, and such a requirement is therefore not needed, so that both cases (absorption and reflection) will be 
treated. Define the constants 


(2.4) 

and the function p,, : ( 0 ,oo) - 


1 


77 := 


2 (p-l)’ 
( 0 , 00 ) by 

\i-p\et """Pi 


p := log(yo) - 


2(l-p) + y2(l-p) 

2iH{l-pY 


et 


{yoy) 


1-p 


(1 _P)2^27 


where I^ is the modified Bessel function of the first kind of order 77 [H Section 9.6]. The CEV density, Cp{y) ■= 
P(Yt G d 77 )/d 77 , then reads 


(2.5) 


p-viy)^ 

Cp(77) = <| v’M 


y^V2Tri 


exp 


(log(p) - yf 

2^H 


if p G [ 1 / 2 , 1 ) or p < 1 with absorption, 
if p > 1 or p < 1 with reflection, 

ifp= 1 , 


valid for y G (0, 00 ). When p > 1, the density converges to zero around the origin, implying that paths are 
being pushed away from the origin. On the other hand diverges to infinity at the origin when p < 1/2, so 
that the paths have a propensity towards the vicinity of the origin. 
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2.3. The moment generating function approach. In the literature on implied volatility asymptotics, the 
moment generating function of the stock price has proved to be an extremely useful tool to obtain sharp 
estimates. This is obviously the case for the wings of the smile (small and large strikes) via Roger Lee’s formula, 
mentioned in Section B.1.11 but also to describe short-and large-maturity asymptotics, as developed for instance 
in m or [39] , via the use of (a refined version of) the Gartner-Ellis theorem. As shown in Section 12.1.11 the 
moment generating function of a stock price satisfying ( 12 . 11 ) is fully determined by that of the random variable V. 
However, even though the density of the latter is known in closed form (Equation (I2.5|) l. the moment generating 
function is not so for general values of p. In the cases p = 0 (with either reflecting or absorbing boundary) and 
p = 1 / 2 , a closed-form expression is available and direct computations are possible. 


2.3.1. Computation of the moment generating function. Denote by Aq Aq ^ and ^12 the moment generating 
function of the random variable V when p = 0 (the subscript ‘r’ / ‘a’ denotes the reflecting / absorbing behaviour 
at the origin) and p = 1/2. Straightforward computations yield 


( 2 . 6 ) 


Kain) 

Ao,r(u) 


= log 


= log 


mt 

I 


2 exp 


{uff^t - 2 ?/o)'« 




_2 
2pou 


exp 


{u^H-2yo)u \ r 2«s/o£; ( + 


+ yo 

yo 


+ -i-s 


V 


fV2t 


- yo 

i^t 


2 - u(H ’ 


where £{z) = exp(—x^)dx is the error function. Note that when p = 1/2 and p = 0 in the absorption 

case, one needs to take into account the mass at zero in (12.31) when computing these expectations. 


2.3.2. Roger Lee’s wing formula. In [45] , Roger Lee provided a precise link between the slope of the total implied 
variance in the wings and the boundaries of the domain of the moment generating function of the stock price. 
More precisely, for any r > 0, let ^^.(t) and U-{t) be defined as 

u+(t) := sup{u > 1 : |A^(u,r)| < 00} and U-{t) ■= sup{u > 0 : |A^(— m,t)| < 00}. 

The implied volatility (Tr(fe) then satisfies 

limsup ^ = '0 (m+(t) — 1) =: /3+(r) and limsup —='0(u_(r))) =:/3_(r), 

fctoo h ki-00 \h\ 

where the function ijj is defined hy if(u) = 2 — 4 — u^. Combining (12. 6 |) and (12.21) yields a closed- 

form expression for the moment generating function of the stock price when p £ {0,1/2}. It is clear that, when 
p = 0 , u±{t) = ±00 for any r > 0 , and hence the slopes of the left and right wings are equal to zero (the 
total variance flattens for small and large strikes). In the case where p = 1/2, explosion will occur as soon as 
{^u{u — l)r^^t — 2 ) = 0 , so that 

u±{t) = i ± and P-{t) = I 3 +{t) = ^ , for all t > 0. 

The left and right slopes are the same, but the product f'^t can be directly calibrated on the observed wings. 
Note that the map r 1 —>■ /3±(t) is concave and increasing from 0 to 2. In [20l [21], the authors highlighted 
some symmetry properties between the small-time behaviour of the smile and its tail asymptotics. We obtain 
here some interesting asymmetry, in the sense that one can observe the same type of rate of explosion (power 
behaviour, given by (ED in the case p < 1), but different tail behaviour for fixed maturity. As r tends to 
infinity, /3 ±(t) converges to 2, so that the implied volatility smile does not ‘flatten out’, as is usually the case 
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for Ito diffusions or affine stochastic volatility models (see for instance [37]). In Section [2751 below, we make 
this more precise by investigating the large-time behaviour of the implied volatility using the density of the 
CEV-distributed variance. 


2.3.3. Small-time asymptotics. In order to study the small-maturity behaviour of the implied volatility, one 
could, whenever the moment generating function of the stock price is available in closed form (e.g. in the case 
p G {0,1/2}), apply the methodology developed in [53]. The latter is based on the Gartner-Ellis theorem, which, 
essentially, consists of finding a smooth convex pointwise limit (as r tends to zero) of some rescaled version of 
the cumulant generating function. In the case where p = 1/2, it is easy to show that 


(2.7) 


lim log 

riO 



- 1 - 00 , otherwise. 


The nature of this limiting behaviour falls outside the scope of the Gartner-Ellis theorem, which requires strict 
convexity of the limiting rescaled cumulant generating function. It is easy to see that any other rescaling would 
yield even more degenerate behaviour. One could adapt the proof of the Gartner-Ellis theorem, as was done 
in [IT] for the small-maturity behaviour of the forward implied volatility smile in the Heston model (see also [T7] 
and references therein for more examples of this kind). In the case (12.71) . we are exactly as in the framework 
of [41) . in which the small-maturity smile (squared) indeed explodes as precisely the same explosion as 

the one in dm. Unfortunately, as we mentioned above, the moment generating function of the stock price is 
not available in general, and this approach is hence not amenable here. 


2.3.4. Large-time asymptotics. The analysis above, based on the moment generating function of the stock price, 
can be carried over to study the large-time behaviour of the latter. In the case p = 1/2, computations are fully 
explicit, and the following pointwise limit follows from simple straightforward manipulations: 

V -l^Z( N / 0: if w 

hm T A (u, r) = < 

Uoo I -|-oo, otherwise. 

The nature of this asymptotic behaviour, again, falls outside the scope of standard large deviations analysis, 
and tedious work, in the spirit of mm, would be needed to pursue this route. 

2.4. Small-time behaviour of option prices and implied volatility. In the Black-Scholes model dSt = 
,/wStdWt starting at So = 1, a European call option with strike e^ and maturity T > 0 is worth 

k \/wT \ kurf ^ VwT^ 


( 2 . 8 ) 


BS(A:,u;,T) = Af - 


\/wT 


-\- 




- eW - 


\/wT 


For any fc S R \ {0}, T > 0, and p > 1, the quantity 


(2.9) 


jP(fc) := < 


BS 


(fc,|,r) y-Pdy, 


if fc > 0, 


^ ^°°(e"-l + Bs(fc,|,r))p-Pdp, ifA:<0, 


is clearly independent of T and is well defined. Indeed, since the stock price is a martingale starting at one. Gall 
options are always bounded above by one, and hence, for fc > 0, J^(fc) < BS{k,y/T,T)y~Pdy -\- y~^dy. 
The second integral is finite since p > 1. When k > 0, the asymptotic behaviour 


BS 


(Lfr) 




,,3/2 
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holds as y tends to zero, so that limj,j,o BS(fc, y/T, T)y ^ = 0, and hence the integral is finite. A similar analysis 
holds when fc < 0 and using put-call parity. Define now the following constants: 


( 2 . 10 ) 


1 


3-2p’ 


Vp ■= 


— p) 


y ■= 


2 


the first two only when p < 1, and note that /3p £ (0,1); define further the following functions from (0, oo) to M: 


( 2 . 11 ) 


/ 

/o(y) 

9o{y) 


e y2(i-p) 

2et(l-p)2’ 

log(y) 

2y et ’ 


/i(y) 

9i{y) 


{yyo)^^ 

log(y) 

et ’ 


as well as the following ones, parameterised by p: 



The following theorem (proved in Section 13.11) is the central result of this paper (although its equivalent 
below, in terms of implied volatility, is more informative for practical purposes): 

Theorem 2.1. The following expansion holds for all k G \ {0} as t tends to zero: 

E (e^- - e'^)’^ = (1 - 6^")+ -f exp - ci(t,p)hi{T,p) + C2(t,p)/i2(T,p)^T'==‘(‘’P^| log(T)|°^(*’^’)c5(t,p) [1 -f 7^(T,p)]. 

Remark 2.2. 

(i) Whenever p < 1, ci and C2 are strictly positive; the function cs is always strictly positive; when p < 1, C3 
is strictly positive; when p = 1, the functions C3 and C4 can take positive and negative values; 

(ii) Whenever p < 1, /i 2 (t, p) < /ii(t, p) for r small enough, so that the leading order is provided by hi] 

(iii) In the lognormal case p = 1, hi{T, 1) ^ (logr)^ as r tends to zero, so that the exponential decay of option 
prices is governed at leading order by exp(—ci(t, l)(logr)^). 

Using Theorem 12.11 and small-maturity asymptotics for the Black-Scholes model (see [551 Corollary 3.5] 
or m). it is straightforward to translate option price asymptotics into asymptotics of the implied volatility: 
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Theorem 2.3. For any fc G K \ {0}, the small-maturity implied volatility smile behaves as follows: 

' k'^^'^t{l — p)' 


al{k) 


(1 - Pp) 
k'^^^t 

rlog(T)2’ 


2 t 


ifP < 1, 
ifp = 1, 

ifp > 1. 


2(2p-l)r|log(r)|’ 

This theorem only presents the leading-order asymptotic behaviour of the implied volatility as the maturity 
becomes small. One could in principle (following [1^ or m) derive higher-order terms, but these additional 
computations would impact the clarity of this singular behaviour. In the at-the-money fc = 0 case, the implied 
volatility converges to a constant: 

Lemma 2.4. The at-the-money implied volatility crT(O) converges to E(-\/V) as t tends to zero. 


The proof of the lemma follows steps analogous to EH Lemma 4.3], and we omit the details here. It in fact 
does not depend on any particular form of distribution of \/V, as long as the expectation exists. Note that, 
from Theorem 12.31 as p approaches I from below, the rate of explosion approaches . When p tends to 1 
from above, the explosion rate is 1 /(t| logrj) instead. So there is a ’’discontinuity” at p = I and the actual rate 
of explosion is less than both these limits. As an immediate consequence of Theorem 12.ll we have the following 
corollary. Define the following functions: 

ci{t,p), ifp<l, 

2 p-I, ifp>I, 

where Ci{t,p) is defined in Table [I] and depends on k (through 



ifp < 1, 

f 

|log(r)-i|, 

if p > 1, 

and A*(fc) := < 

log(T)-^ 

ifp= 1 , 

1 


Corollary 2.5. For any p £ K, the sequence satisfies a large deviations principle with speed h*{T,p) 

and rate function A* as r tends to zero. Furthermore, the rate function is good only when p < 1. 

Recall that a real-valued sequence (3n)n>o satisfies a large deviations principle (see [18] for a precise intro¬ 
duction to the topic) with speed n and rate function A* if, for any Borel subset B C U, the inequalities 

— inf A*(z) < liminf logP(3n G B) < limsupn”^ logP(3n G B) < — inHA*(2:) 

2 GB° nfoc nfoc zGB 

hold, where B and B° denote the closure and interior of B in K. The rate function A* : K. —?► K U {-boo}, by 
definition, is a lower semi-continuous, non-negative and not identically infinite, function such that the level sets 
{cc G M : A*(a:) < a} are closed for all a > 0. It is said to be a good rate function when these level sets are 
compact (in R). 


Proof. The proof of Theorem 12.11 holds with only minor modifications for digital options, which are equivalent 
to probabilities of the form P {Zr < fc) or P {Zr > fc). For p £ (—oo, 1], one can then show that 

hmfc*(T,p) logF {Zr < k) = — inf {A*(x) : x < fc} . 

The infimum is null whenever fc > 0 and p < 1, and A}(x) = l/(2^^t) is constant. Consider now an open 
interval (a, b) C R. Since (a, b) = (—oo, b) \ (—oo, a], then by continuity and convexity of A*, we obtain 

lim /i*(t,p) logP (Z,- £ (a, 6)) = — inf A*(x). 

riO xG{a,b) ‘ 
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Since any Borel set of the real line can be written as a (countable) union / intersection of open intervals, the 
corollary follows from the definition of the large deviations principle [HI Section 1.2]. Whenp S (l,oo), the only 
non-trivial choice of speed is |(logr)“^|, in which case limT^o l(logT)“^| \og¥ {Zr < k) = —{2p— 1). Clearly, the 
constant function is a rate function (the level sets, either the empty set or the real line, being closed in K), and 
the corollary follows. □ 

Remark 2.6. In the case p = 1/2, as discussed in Section [2.3.4l the moment generating function of Z is available 
in closed form. However, the large deviations principle does not follow from the Gartner-Ellis theorem, since 
the pointwise rescaled limit of the mgf is degenerate (in the sense of (12.71) 1. 


2.4.1. Small-maturity at-the-money skew and convexity. The goal of this section is to compute asymptotics 
for the at-the-money skew and convexity of the smile as the maturity becomes small. These quantities are 
useful to traders who actually observe them (or approximations thereof) on real data. We define the left 
and right derivatives by i9/"cr^(0) := limfc-|-o 9fccr^(fc)|fe=o and i9/'fT^(0) := limfe|o'9feO'?(^)|fc=o, and similarly 
^fefeO-?(0) := limfcto 5fcfcO'r(fc)|fe=o and cI/).(t^( 0) := limfe^o i9fefe(T^(fc)|fe=o- The following lemma describes this 
short-maturity behaviour in the general case where V is any random variable supported on [0,oo). 

Lemma 2.7. Consider (12.11) and assume that E(V"/^) < oo for n = —1,1,3, and uit := P(V = 0) < 1. As t 
tends to zero, 

dkk^riO) 

Jensen’s inequality and the fact that the support of V is in [0,oo) imply that both E(V^/^) — E(-\/V)^ 
and E(V-i/2) -E(v^)-i are strictly positive. The small-maturity at-the-money skew is always negative for 
small mt. This in particular means that the smile generated by m is not necessarily symmetric. When m* > 0, 
the at-the-money left skew explodes to —oo and the at-the-money right skew explodes to -foo. Furthermore, 
the small-maturity at-the-money convexity tends to infinity. 


Kf’Uo) ~ s/s (e (v-‘/“) - e(vT)-' . 


Proof. We first focus on the at-the-money skew. By definition C{k,T) = BS(fc, CT^(fc), r) and therefore 

dkC{k,T) = 5feBS(A:,cr/(A:),r) -H dkcrl{k)d^uBS{k,a'^{k),T). 

Also by (13.11) , an immediate application of Fubini yields 

pOO 

(2.12) dkC{k,r)= / dkBS{k,y,r)Cp{y)dy-\-mtdk{l-e^)~^, 

Jo 

We first assume that mt = 0. The at-the-money skew is then given by 

9fccr^(fc)|fe=o = {dn,BS{k,al{k),T)\k=o) ^ (^J 9fcBS(fc,p,r)|fc=oCp(y)dy - 9fcBS(fc,cr^(fc),r)|fe=o^ ■ 

Recall now from Lemma [T4l that 17^(0) = E(-\/V) -I- o(l). Straightforward computations then yield 

T3('n^^-3/2 


(2.13) 


^ -V (^) . 4 + . o ( 4 ( 0 ,.-) , 


5„BS(fc,<(fc),T)|fe=o = 


tTr 


(O)v^ 


exp 
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as T tends to zero. Hence 

dk(rl{k)\k=o = exp ^]E(’/V) - crT(O) 

and so, as T tends to zero, 


(E(V3/2) -cr^(0)3)7 
24 


OK(0) 


V)). 


(2.14) ~ -5^ (e(V’/“) - (E(A1))“) t, 

The small-maturity convexity follows similar arguments, which we only outline: 

(2.15) dkkC{k,T) = i9fefeBS(fc, tT^(fc),r) -h 2afc(T^(A:)5u,feBS(fc, cr^(fc), r) -h {dkcrl{k)Y i9u,u,BS(fc, CT^(fc), r) 

dkk(jl{k)dwBS{k, cjl{k), r). 


and 

(2.16) 


poo 

dkkC{k,r)= / «9fcfcBS(A:,y,r)Cp(y)dy + mf5fcfc (l-e^) 
Jo 


Likewise, we first consider the case mj = 0. Straightforward computations yield 


(2.17) 


exp ( - 

afcfcBS(fc,tT^(fc),T)|fc=0 = 

9feu,BS(fc,cr^(fc),r)|fc=o = exp ( - 

a^„BS(fc,cr^(fc),T)|fc=o = - 


(_iMi) 


(Tr{0)y/T\/^ 


-M 


^ crr(0)v^ ^ ^ 


r^e 


8 / 4cri-(0)-\/^’ 

'l^(ra2(0)+4) 


v^ctt-(0 )i/r 2 

y/r Sr^/^ariO) 


- o +c> K (0)4/7), 


o 


7-3/2 

’■r(O) 


16-\/27(t(t^( 0))^/^ 4-\/27rtT3 (0) 512-\/27 

Using (I2.15I) and (I2.16I) in conjunction with (I2.13l) . (l2.l7l) and ()2.14p . we obtain dhha-(O) ^ i]E(-\/V) — E(-\/V)“^y 

When mt > 0, we need to take right and left derivatives in (I2.12I) and (I2.16I) to account for the atomic term. 

Since (l - e'^)"'’ |fc=o = d^k lfc=o = “1 and 9/ (l - e'')'*’ |fc=o = d^k (l “ lfc=o = 0, the lemma 

follows immediately. □ 


2.5. Large-time behaviour of option prices and implied volatility. In this section we compute the large¬ 
time behaviour of option prices and implied volatility. The proofs are given in Section 13.21 It turns out that 
asymptotics are degenerate in the sense that option prices decay algebraically to their intrinsic values. The 
structure of the asymptotic depends on the CEV parameter p and whether the origin is reflecting or absorbing: 


Theorem 2.8. Define the following quantity: 

«„) .=_ 24-«^-4r(i-2p) _ / vl"--' \ 

with r] given in The following expansions hold for all k € as t tends to infinity: 

(i) if p <3/4: and the origin is absorbing then 

E (e^- - e'=)+ = 1 - mt + mt(l - e^)+ - Se'^/^yo Q - 

(ii) if p <1/4 and the origin is reflecting then 

E (e^- - e'=)+ = 1 - ^ XL ■ 
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For other values of p, asymptotics are more difficult to derive and we leave this for future research. The 
asymptotic behaviour of option prices is fundamentally different to Black-Scholes asymptotics (Lemma IA.2|) 
and it is not clear that one can deduce asymptotics for the implied volatility. For example, the intrinsic values 
do not necessarily match as r tends to infinity because of the mass at the origin. The one exception is when 
the origin is reflecting, in which case the implied volatility tends to zero. This result follows directly from the 
comparison of Theorem 12.81 and Lemma IA.2I 

Theorem 2.9. If p <1/4 and the origin is reflecting, then for all k G R as t tends to infinity: 

aUk) ~ 

Although, we have provided the large-time asymptotics in this section, it is not our intention to use this model 
for options with large expiries. Our intention (as mentioned in Section[T]) is to use these models as building blocks 
for more complicated models (such as stochastic volatility models where the initial variance is sampled from a 
continuous distribution) so that we are able to better match steep small-maturity observed smiles. In these types 
of more sophisticated models, the large-time behaviour is governed more from the chosen stochastic volatility 
model rather than the choice of distribution for the initial variance (see [3H1 HO] for examples), especially if the 
variance process possesses some ergodic properties. This also suggests to use this class of models to introduce 
two different time scales: one to match the small-time smile (the distribution for the initial variance) and one 
to match the medium- to large-time smile (the chosen stochastic volatility model). 

2.6. Numerics. We calculate option prices using the representation (EH) and a global adaptive Gauss-Kronrod 
quadrature scheme. We then compute the smile with a simple root-finding algorithm. In Figure [T] we plot 
the smile for different maturities and values for the CEV power p. The model parameters are yo = 0.07, 
f = 0.2yy^~^ and t = 1/2. Note here that we set ^ to be a different value for each p. This is done so that 
the models are comparable: f is then given in the same units and the quadratic variation of the CEV variance 
dynamics are approximately matched for different values of p. The graphs highlight the steepness of the smiles 
as the maturity gets smaller and the role of p in the shape of the small-maturity smile. Out-of-the money 
volatilities (for K ^ [0.9,1.1]) explode at a quicker rate as p increases (this can be seen from Theorem l2.3|) . The 
volatility for strikes close to the money K £ [0.9,1.1] appears to be less explosive as one increases p. This can 
be explained from the strike dependence of the coefficients of the asymptotic in Theorem 12.31 and is discussed 
further in Section [577 t1 

2.7. Application to forward smile asymptotics. We now show how our model EH and the asymptotics 
derived above for the implied volatility can be directly translated into asymptotics of the forward implied 
volatility in stochastic volatility models. For a given martingale process e^, a forward-start option with reset 
date t, maturity t + t and strike e^' is worth, at inception, E(exp(At+T- — Xt) — e^)+. In the Black-Scholes 
model, the stationarity property of the increments imply that this option is simply equal to a standard Call 
option on e'’^ (started at Xq = 0) with strike e^ and maturity r; therefore, one can define the forward implied 
volatility crt,T{k), similarly to the standard implied volatility (see |39| for more details). Suppose now that the 
log stock price process X satisfies the following SDE: 

dXs = ——Ygds + Aq = 0 , 
di; = fsYfdBs, Yo = yo> 0, 

d(W^,S)^ =/3ds. 


( 2 . 18 ) 
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Figure 1. Here we plot K i—>■ crr{\ogK) for maturities of 1 (circles), 1/2 (squares) and 1/12 
(diamonds) for increasing values of the CEV power p. The smile is obtained by numerically 
solving for the option price using (13.111 and then using a simple root search to solve for the 
implied volatility. Parameters of the model are given in the text. 


with p £ R, IpI < 1 and W, B are two standard Brownian motions. Fix the forward-start date t > 0 and set 


(2.19) 


in 


if 0 < M < t, 
if u > t, 


where ^ > 0 and ^ > 0. This includes the Heston model and 3/2 model with zero mean reversion (p = 1/2 
and p = 3/2 respectively) as well as the SABR model (p = 1). Consider the CEV process for the variance: 
dVu = To = 2/0, where p £ M and H is a standard Brownian motion. Let Xr^ := Xt+r — Xr denote 

the forward price process and let CEV(t,^,p) be the distribution such that Law(Yt) = Law(V) = CEV(t,^,p). 
Then the following lemma holds: 


Lemma 2.10. In the model (I2.18P the forward price process 


solves the following system of SDEs: 


( 2 . 20 ) 


dXi*'> = -f \/^dWr, = 0, 

dV-(‘) = e dBr, - CEV(C e, p), 

d (VF, B)^ = pdr, 


where Yr 


(t) 


is independent to the Brownian motions (lFT)r>o (i3r)r>o- 
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This lemma makes it clear that forward-start options in stochastic volatility models are European options on 
a stock price with similar dynamics to (12.181) . but with initial variance sampled from the variance distribution 
at the forward-start date. If we set ^ = 0, then = Z and the following corollary provides forward smile 
asymptotics: 

Corollary 2.11. If^=0, Theorem \2.1[ Theorem \2.S\ and Lemma \2.4\ hold with Z = and CTt = cr-t T. 

Remark 2.12. 

(i) This result explicitly links the shape and fatness of the right tail of the variance distribution at the 
forward-start date and the asymptotic form and explosion rate of the small-maturity forward smile. Take 
for example p > 1: the density of the variance in the right wing is dominated by the polynomial 

and the exponential dependence on y is irrelevant. So the smaller p in this case, the fatter the right tail 
and hence the larger the coefficient of the expansion. This also explains the algebraic (not exponential) 
dependence for forward-start option prices. 

(ii) The asymptotics in the p > 1 case are extreme and the algebraic dependence on r is similar to small- 
maturity exponential Levy models. This extreme nature is related to the fatness of the right tail of the 
variance distribution: for example, the 3/2 model {p = 3/2) allows for the occurrence of extreme paths 
with periods of very high instantaneous volatility (see [HJ Figure 3 ]). 

(iii) The asymptotics in Theorems l2.1l and l2.3l remain the same (at this order) regardless of whether the variance 
process is absorbing or reflecting at zero when p G (—oo,l/2). Intuitively this is because absorption or 
reflection primarily influences the left tail whereas small-maturity forward smile asymptotics are influenced 
by the shape of the right tail of the variance distribution. 

2.7.1. Conjecture. Whenp = 1/2 in Corollarv l2.111 the asymptotics are the same as in [HI Theorem 4.1] for the 
Heston model. This confirms that the key quantity determining the small-maturity forward smile explosion rate 
is the variance distribution at the forward-start date. The dynamics of the stock price are actually irrelevant 
at this order. This leads us to the following conjecture: 

Conjecture 2.13. The leading-order small-maturity forward smile asymptotics generated from (I2.18P are equiv¬ 
alent to those given in Corollaru \2.11[ 

Practitioners have stated nng that the Heston model (p = 1/2) produces small-maturity forward smiles 
that are too convex and ” U-shaped” and inconsistent with observations. Furthermore, it has been empirically 
stated m that SABR or lognormal based models for the variance (p = 1) produce less convex or ” U-shaped” 
small-maturity forward smiles. Our results provide theoretical insight into this effect. We observed in Section|2]6| 
and Figure |T] that the explosion effect was more stable for strikes close to the money as one increased p. The 
strike dependence of the asymptotic implied volatility in Theorem 12.31 is given by iL i—>■ a/| logiL| for p = 1/2 
and K H> | logitr| for p = 1. It is clear that the shape of the forward implied volatility is more stable and less 
U-shaped in the lognormal p = 1 case. 

3. Proofs 

3.1. Proof of Theorem 12.11 Let C{k, r) := E(e^^ — e^')+. This function clearly depends on the parameter t, 
but we omit this dependence in the notations. The tower property implies 

BS{k,y,T)(p{y)dy -\-mt (l-e'^)^. 



(3.1) 


C{k,r) 











BLACK-SCHOLES IN A CEV RANDOM ENVIRONMENT: A NEW APPROACH TO SMILE MODELLING 


15 


where BS is defined in (12.81) . Cp is density of V given in (12.5|) and m* is the mass at the origin (12.31) . Our goal 
is to understand the asymptotics of this integral as r tends to zero. We break the proof of Theorem 12.11 into 
three parts: in Section 13.1.11 we prove the case p > 1, in Section 13.1.21 we prove the case p G (—oo, 1) and in 
Section 13.1.31 we prove the case p = 1. We only prove the result for fc > 0, the arguments being completely 
analogous when fc < 0. The key insight is that one has to re-scale the variance in terms of the maturity r 
before asymptotics can be computed. The nature of the re-scaling depends critically on the CEV power p and 
fundamentally different asymptotics result in each case. Note that for fc > 0, (l — e^)^ = 0, so that the atomic 
term in dSU is irrelevant for the analysis. When fc < 0 the arguments are analogous by Put-Call parity. 


3.1.1. Case: p > 1. In Lemma [XT] we prove a bound on the CEV density. This is sufficient to allow us to prove 
asymptotics for option prices in Lemma 13.21 after rescaling the variance by r. This rescaling is critical because 
it is the only one making JiS{k,y /t,t) independent of r. Let 


X{r,p) ■■= --- , 

|1 -p|^2^r(l -b |77|)(^2(1 -p)2^2i 

and we have the following lemma: 


kl 


exp - 




Lemma 3.1. The following bounds hold for the CEV density for all y,T > 0 when p > 1: 


(-) 

^ xicp) j 

\r/ 

y^^ \ 


/ x(cp) j 

Vr / 

- y 2 p 1 


1 - 


1 


2et{l-pr \y 

1 + exp 


2p-2' 


2-2p 

^0 


2{p-me 


2et{i-pr \y 


2p-2 


^2t(l -_p)2 (yyo) 


P-1 


Proof. From ill Equation (6.25)] we know that for a; > 0 and v > —1/2, the modified Bessel function satisfies 


(3.2) 


r(i. + i) 


/x\'' fx\ 

(2) s ^ i>TI) (5) 


so that the expression for the CEV density in (12.51) implies that for p > 1, 


Xij,p) ( 1 f t\ 

2et{i-pr[y) 


r/2p 


2p-2' 


< 




where 


^{y,T) := - 


1 


2f2t(i-_p)2 


2p-2 


y2p 

r 

yyo 


p-i 


^2t(l -p)2 

For fixed r > 0, note that m{-,T) : M+ 1—> K+ takes a maximum positive value at p = y^r with m{yoT,T) = 
Pq~^^/(2(p — l)'^tf'^). When m(-) > 0 Taylor’s Theorem with remainder yields = 1 -b e^m{y,T) for 

some 7 G (0,m(p,T)), and hence < 1 -b e^^y°'^’'^'^m{y,T). If m(-) < 0 then < I -b |m(p,T)| < 

1 -b e'"*^^“'’’’'’’^|TO(p, t)|. The result for the upper bound then follows by the triangle inequality for |to(p, t)|. The 
lower bound simply follows from the inequality I — x < , valid for x > 0, and 


1 - 


2fH{l-pY \y 


2p-2 


< exp — 


2f‘^t{l -p)2 \y 


2p-2^ 


□ 


Lemma 3.2. When p > Theorem \2.1\ holds. 
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Proof. The substitution y ^ y/r and (13.11) imply that the option price reads C(fc, r) = BS(fc, y, T)fp{y)dy = 
T~^ BS{k,y/T,T)(p{y/T)dy. Using Lemma [XT] and Definition (12.91) . we obtain the following bounds: 


Xir,p) 


3^P{k) - 


^ 2 p -2 


Xir,p) 


j 2 P(fc) + exp 


2-2p 

Vo 


r 2 p -2 


.J4p-2(fc) + 


2fH{l-p)' 

T-P-l 


:3^P-^{k) 


-3^P-\k) 


Hence for r < 1: 

C{k, t)t 


x(r,p)J 2 p(fc) 


- 1 


2{p-l)He) \2et{l-pY'- ' ^2t(i-p)2yP 

( j4P-2(/c) j3P-l(fc) 


<C'(fc,T), 

> C'(fc,r). 


< exp 


2-2p 

2/0 


2{p - l)He I \ 2et{l - p)2 J2p(fc) ^2^(1 _ p)2yP-l J2p(fc) 


rP -1 


which proves the lemma since 3'^P{k) is strictly positive, finite and independent of r. 


□ 


3.1.2. Case: p < 1. We use the representation in (13.11) and break the domain of the integral up into a compact 
part and an infinite (tail) one. We prove in Lemma 13.41 that the tail integral is exponentially sub-dominant 
(compared to the compact part) and derive asymptotics for the integral in Lemma [331 This allows us to apply 
the Laplace method to the integral. We start with the following bound for the modified Bessel function of the 
first kind and then prove a tail estimate in Lemma 13.41 


Lemma 3.3. The following bound holds for all x > 0 and v > —3/2; 

(I)' 


T ; \ ^ + 2 (X\- 2x 

luix) < ——(-1 


T{u + 2) 

Proof. Let x > 0. From [55l Theorem 7, page 522], the inequality lu{x) < fo+i(x)^/U+ 2 (x) holds whenever 
12 > —2, and hence combining it with (13.21) (valid only for v > —1/2), we can write 

r(i/-|-3) /xN^ 


(f) 


when V > —3/2. The lemma then follows from the trivial identity Tiy -|- 3) = (i 2 + 2 )r(i 2 -|- 2). 
Lemma 3.4. Let L > 1 and p < 1. Then the following tail estimate holds as t tends to zero: 


□ 


- p) 


I Bs(fc,^,r)Cp(^)d2/ = o(exp(-^ 
Proof. Lemma 13.31 and the density in (12.51) imply 


T 

^ 1-p 

4l-/3p)/2 yo 


2\ 


^0 _ 2 p / 1 \ ^ 

:^y 


2f,'^t{l — p)"^ 


1—p 

Vo 


( 2 / 2 / 0 ) 


i-p 


r^p(i Pl^‘^t{l — p)'^ 


where the constant bo is given by 

(17 + 2 ) 


|i-p|^2ir(^ + 2)(2(i-p)2e2ty 


resp. 


(H + 2) 


II- 


p\etT{\y\ + 2)(2{l-p)^et 


\v\' 
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if the origin is reflecting (resp. absorbing) when p < 1 / 2 ; the exact value of bo is however irrelevant for the 
analysis. Set now L > 1. Using this upper bound and the no-arbitrage inequality BS(-) < 1, we find 


BS 

bo 


{k.^,r)c,{:^)dy< [ Cp(^)dy 


T> 

/*oo 


< 

^—2p^p 


y ^^exp 


y^ ^^exp 


1 


— p)'^ 
1 


i-p 

- Vo 


(yya) 


i-p 


r^p(i P)^'^t{l—p)^ 


2^H{l-py \r/3p(i-p) 


y^-^ i-P 

- yo 


(.yyo) 


i-p 


r/5p(i p)^‘^t(l — p)'^ 


dy 


dy, 


where the last line follows since y^ > y Setting q = ^y^ Pjrl^pi'^ p) — p^ /(^\/t(l — p)) yields 

dy 



/ 

( ^ 


2 \ 

y^ exp 



yo ^ 

(yyo)^ ^ ] 

1 

2 ^ 2 i(l _p )2 

r^p(i“P)^ 2 i(x — p )2 1 


(3.3) = 


^Vt(i-p) 

'j-2j3p {p— 1) 


pOO 

^Vi{l-p) / q exp 

JLr 


2 1 —P 

9 , yo 9 


2 ^v^(l-p) 


dq + yl ^ I exp 


pOO 

Jl^ 


2 1 —P 

9 , % 9 


2 ^Vi(l-p) 


dq 


with Lr ■= Pjrl^pi'^ p) — p^ /(^y/i(l — p)) >0 for small enough r since L > 1 and p G (—oo, 1). Set 
now (we always choose the positive root) 


L^-P 


(/3p(l-p)) 


so 


that, for T < T* we have Lr > 4yl ^/(^y/i{l — p)) and hence for q > Lr'. 


1—p 2 

Vo 9 ^ 9 


^Vt(l-p 4 

In particular, for the integrals in (13.31) we have the following bounds for t < t*: 


L 


2 1—P 

9 , 90 9 


2 ^v^(l-p) 


dq< q exp —— ) dg = 2 exp — 


2 1 —P 

9 , 9o ^9 


2 ^^/i(l-p) 


dq < J exp ( —y ) dg < -^ exp ( —j- ) : 


where the last inequality follows from the upper bound for the complementary normal distribution function 
in [591 Section 14.8]. The lemma then follows from noting that 1 — /3p = 2/3p(l — p). □ 

Lemma 3.5. When p < 1, Theorem \2.1\ holds. 

Proof. Let f := with (3p defined in (12.101) . Applying the substitution j/ —>• y/r to (13.11) yields 
C{k,T)=J BS(fc,y,T)Cp(?/)dy = BS |,r) Cp (|) dy 




BS 


(/c, |, r) Cp (I) dy + ^ y BS (fc, |, r) Cp (|) dy. 
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for some L > 0 to be chosen later. We start with the first integral. Using the asymptotics for the modified 
Bessel function of the first kind [I] Section 9.7.1] as r tends to zero, we obtain 


m) 


(-3p^p/2yP/2g~ 252 t(i_p)i! 


^y3p/2.^27rt 


■ exp - 


1 


y 


2(l-p) 


1 {yyoY^ 

r2/3p(i-p) 2^H{1 - p)2 '■ r/5p(i-p) ^ 2^(1 _ p )2 


1 + O 


Note that this expansion does not depend on the sign of rj and so the same asymptotics hold regardless of whether 
the origin is reflecting or absorbing. In the Black-Scholes model, Call option prices satisfy fLemma lA.ll) : 


BS 


(fc,|,r) 


,,3/2 /-rx3/2 


exp 


T k 


as T tends to zero. Using the identity 1 — /3p = 2/3p(l — p) we then compute 






_r'o_L 




% e 


o . , /o(i/) 1 /i(y) 

/ yifi-P/e 

1 + 0 yp)/2'\ 

Jo 

L V /J 


2t: k'^^y/t 

where /o, /i are dehned in (12.111) . Solving the equation Jq^u) = 0 gives y = yp with Pp defined in (12.101) and we 
always choose the positive root and set L > Pp. Let /(t) := exp + ^(/.igp)/2 ^ dy. Then for 

some e > 0 small enough, as t tends to zero: 


/(r)-exp --j^ + 


foiPp) , fiiPp) , fiiyp)^ 


t1 r(i-^p)/2 2f'^{y ) 


yp 


/•Pp + E 

/ 

1 

/ exp 


2 

Jpp-e 

V 



fo{yp){y-yp) f[{y\ 


'-(1 ^p)/2 


/o'(2/p) 


dy 


+ 2f''{yp) J 

12 ' 


exp 


V 


/o'(2/p)(2/-2/p) /](m) 


n(l—/3p)/2 


/^(yp) 


dy 


= exp 


fo(yp) fijyp) fiiyp? \ _a-3.)/2-j(i-p) / 2^ 


T-l /3p .p(l-/3p)/2 2fll{y ) 


y0 


fHiypY 


The ~ approximations here are exactly of the same type as in |32j . and we refer the interested reader to this 
paper for details. It follows that as r tends to zero: 

7^: / BS (k, r) Cp dy = exp C5(t,p)T-(‘-rf [l + O (r^) 


From Lemma 13.41 we know that 


y BS {k, ■^,T^ Cp(y//3p)dy = O f exp 


L^-P 


2et{i-p) \Ti^-Pp)n 


% ^ 


Choosing L > max ^1, (2f^t(l — p)/o(yp)) ^^\yp'^ makes this tail term exponentially subdominant to 

/o^ BS(fc,y/T^p,T)Cp(2///3p)dy, which completes the proof of the lemma. □ 
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3.1.3. Case: p = \. We now consider the lognormal case p = 1. In this case, the random variable log(V) is 
Gaussian with mean p (defined in (I2.4|) l and variance The proof is similar to Section [3.1.21 but we need 
to re-scale the variance by t| log(T)|. We prove a tail estimate in Lemma 13.61 and derive asymptotics for option 
prices in Lemma 13.71 


Lemma 3.6. The following tail estimate holds for p = 1 and L > 0 as t tends to zero (p defined in (12.4p ).’ 

n 2l \ 


/: 


BS k, 


y 


r|log(T)|^ 


Cl 


y 


T|l0g(T)| 


dy = O \ exp < —; 


1 


log 


L 


L V'r|log(r)| 




Proof. By no-arbitrage arguments, the Call price is always bounded above by one, so that 


BS k 


y 


’ Tiiog(T)r 


Cl 


y 


T|log(T)| 


dy < Cl 


y 


r|log(T)| 


dy. 


With the substitution q = ^^[log(y/(r| log(r)|)) —y], the lemma follows from the bound for the complementary 
Gaussian distribution function [59l Section 14.8]. □ 

Lemma 3.7. Let p = 1. The following expansion holds for option prices as t tends to zero: 

C{k,T) = C 5 {t, l)exp - ci{t, l)hi{T,p) + C2(t, | log(T)|°^(*’^^ ^1 -b O 

with the functions ci, C2,..., C5, hi and h^ given in Table\^ 


log(r)| 


Proof. Let r := rj log(T)|. With the substitution y ^ y/r and using (13.11) . the option price is given by 

pOO -t pOO 

C{k,T) = J BS(fc,y,T)Ci(y)dy = BS |,rj Ci (|) dy 


1 

T 


BS 


(k, t) Cl (I) ^y + J I’ '^) Cl (I) dyj =: C(fc, t) + C{k, r). 


for some L > 0. Consider the first term. Using Lemma [A. II with r = rj log(r)|, we have, as r tends to zero. 


BS k 


’ 'r|l0g(T)p 


, T = exp 


fc2|log(r)| k 


.,3/2 


Therefore 
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I l0g(T)|3/2^fc227ry^ /q 


exp 


2y 2/ fc^l log(r)p/2-y^ 

^ -. (log( " 


1 -bO 
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fe^|log(T)| 

2y 


r| Iog(T)| 


2^H 


log(r)| 


y^/^dy 


fk (log(T)-blog] log(r) 1)2-b/i^ 

/i(log(T) +log|log(r)|)^ 

h{T) 

[l+^(|log(r)|)] 

\2 2Ct 

et y 

^A:227r-\/t| log(r)p/2 


where /i(t) := g2{y) exp (-yo(y)l logrj -b yi(y) log | log(T)|) dy with yo and gi defined in (12.111) and 


52 (y) := Vyexp 


Mlog(y)\ 


et )■ 

The dominant contribution from the integrand is the j log(r)| term; the minimum of yo is attained at y* given 
in (12.101) . and yo(y*) = > 0. Set 

2N 

dy = 


Io{t):=J exp(-i |^(y - y*)Y^| log(T)|y"(y*) - 


y'(y*)log|log(T)| 


27r 


v'|log(T)|y(,'(y* 


y"(y*)|log(r)|- 
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Then for some e > 0 as r tends to zero, with L > y*, 
/•y'+e 


h{T) 


fy*-e 


52(2/) exp I - 30(2/)! log(r)I + 51 ( 2 /) log|log(T)||d2/ 

ry 

^2(2/*)e“®°l^*ll '°s('r)l+9i(y*) '°s I iog(T)| / (y*)(y-y*)^l iog('r)l+9i(y*)(y-y*) log I logi'^il^jy 

Jy*-e 

W,{y*)\og\\og{T)\f 


~ 52(2/*) exp (^-£/o(2/*)|log(r)| + 5i(2/*) log | log(T)| 
= 92 {y*) exp {-go{y *) | log(T) | + gi {y *) log | log(T) | 


2g”{y*)\\og{T)\ 

{g[{y*)\og\\og{r)\)^ 


Io{t) 


27r 


25o ( 2/*)l logMI J ]] goiy*)\ log(r)|' 

where again the ^ approximations here are exactly of the same type as in |32) . and we refer the interested 
reader to this paper for details. Therefore as r tends to zero: 

1 


C(fc, r) = C 5 {t, 1) exp - ci(t, 1 )/ii(t, 1) + C2(t, 1 )/i 2 (t, 1))| log(r)|‘='‘l*’^l 


1 + 0 


|log(T)|yJ ’ 

with the functions ci, C2,C5, hi and /i2 given in Table [TJ For ease of computation we note that 




C5(t, 1) = 

Now by Lemma 13.61 


Miog(y*) 

iH 


) _ l^l^^^^/2 exp(|-, ^ 2 * 


Miog(y*)' 


C(fc,r) = 


|log(T)| 


fc2^v^\/go(2/* 


BS k, 


40r 
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r|log(T)| 
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r|log(r)| 


O exp <^ - 


r|log(T)| 
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log 


L 


r|log(r)| 


1 2 


Since for some B > 0 we have that 
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exp - 


2^H 


log 


L 


I log (t) I 


1 2N 




<B(r|log(T)|)i^ 


ir(log(T)-M) 


exp(-^/ii(Tl) 


choosing L such that log(L) > y yields 

1 


O exp<^ -- 


log 




2^H [ V'l'IlogMI 

Hence C{k,T) is then exponentially subdominant to the compact part since 


exp 1 )/ii(t, 1) - C2(t, 1 )/i2(t, l)jO |^exp ■ 

and the result follows. 


1 


2^H 


‘“Km)-"] 


□ 


3.2. Proof of Theorem 12.81 Lemma [4.21 and (13.11) yield the following asymptotics as r tends to infinity: 
(3.4) C{k, r) = 1 - mt + mt(l - e^')+ + 

where i2(T) := q{z)e~'^^dz, and we set q{z) = —8(^p{8z)/^/¥z. As z tends to zero recall the following 

asymptotics for the modihed Bessel function of the first kind of order 77 [TJ Section 9.6.10]: 
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Using this asymptotic and the definition of the density in (12.51) we obtain the following asymptotics for the 
density as y tends to zero when p < 1 and absorption at the origin when p < 1/2; 


(3.5) 


Cpiy) = 


yoy 


l-2p 


1 - p|c^tr(|?7| +1) (2(1 - py^'^t) 


hi 


exp - 


2et{i-py 




Analogous arguments yield that when p < 1/2 and the origin is reflecting, then, as y tends to zero. 


(3.6) 


Cpiy) 


|i_p|^2ir(p + i)(2(i-p)W 



2^H{i-py 




In order to apply Watson’s lemma m Part 2, Chapter 2] to £, it is sufficient to require that the function q 
satisfies q{z) = 0(e^^) for some c > 0 as z tends to infinity. This clearly holds here since lim^ioo Cp(^) = 0. We 
also require that q(z) = aoz^(l + ©(z”)) as z tends to zero for some I > —1 and n > 0. When p > 1, it can be 
shown that Cp is exponentially small, and a different method needs to be used. When p < 1 and the density is 
as in (13.51) then I = 1 — 2p — 1/2 and so we require p < 3/4. Analogously, when p < 1/2 and the density is (13.6|) 
then I = —2p — 1/2 and we require p < 1/4. An application of Watson’s Lemma in conjunction with (13.41) then 
yields Theorem 12.81 


Appendix A. Black-Scholes asymptotics 

This appendix gathers some useful expansions for the Black-Scholes Call price function BS defined in (12.81) . 

Lemma A.l. Let k,y > 0 and r : (0, oo) —>• (0, oo) be a eontinuous function sueh that lim ^ = 0. Then 

T^o t(t) 


BS k 


,,3/2 


(- 


t(t) ’ J k'^\/^ \TiT) 


•ill 


exp - 


k^ r(r) k 
-LA + _ 

2y T 2 


1 + 0 


^-ir) 


as T tends to zero. 


Proof. Let k,y > 0 and set t*(t) = t/t(t). By assumption, r*(r) tends to zero as t approaches zero, and 
hence (12.81) implies 

BS {k, |,t) = BS {k,y,T*{T)) = ^^{d*_^_{T)) - e''A/'(dl(r)), 

where we set d’^{T) := —k/{y^yT*{T)) ± i\/pT*(r), and Af is the standard normal distribution function. Note 
that tends to —oo as r tends to zero. The asymptotic expansion 1—Af{z) = {2n)~^/^e~^^ (z“^ — z~^ -I- 0(z ~^)), 
valid for large z m page 932]), yields 


BS 


=Af(d;(r)) - e^Af {d*_iT)) = 1 - AT (-d+(T)) - e'=(l - Ar(-d!L(T))) 




exp --dl(T)V2 
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1 


dl(r) d;(r) d;(r)3 dl(r)3 

¥1 


+ 0 


d;(r)5 


as T tends to zero, where we used the identity |(il(T)^ — k = ^d’^iry. The lemma then follows from the 
following expansions as r tends to zero: 
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2-+V 


exp ( -irdiyy ) = exp ( - 
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□ 



























22 


ANTOINE JACQUIER AND PATRICK ROOME 


Lemma A.2. Let y > 0 and /c G R. Then, as t tends to infinity, 

‘ (“f + f) ■ 

Proof. Let y > 0. Then 

BS(fc,y,T) = Af [d*_^_{T)) - e’^Af {d*_{T)) , 

where we set d’^{T) := —k/{y/^) ± ^y/yr^ is the standard normal distribution function. Hence dj. tends 

to ±oo as T tends to infinity. The asymptotic expansion 1 — N{z) = (27r)“^/^e“^^/^ (z“^ — z~^ + 0{z~^)), 
valid for large z m page 932]), yields 

BS(fe, y,T)=Ar «(r)) - e'^ (1 - AT {-d*_ (r))) 


= 1 - 


1 




exp |- 


1 


1 


o 


1 


d;(r)5 


.d;(T) d*_(r) dl(r)3 d;(T)3 

as r tends to infinity, where we used the identity ^d’f{T)^ — k = ^dii_(r)^. The lemma then follows from the 
following expansions as r tends to inhnity: 


exp = exp (i + 0{t i)) , 


d+(T) dl(T) dl(T)3 d+(r)3 y/TFry 


(l + 0(r-i)). 


□ 


Appendix B. Hedging 


We suppose that the market hltration (A't)t>o does not include cr(V); otherwise hedging reduces to the Black- 
Scholes case. We also suppose for simplicity that P(V = 0) = 0. We set Ct{T, k) := £((51- — e^')+|A't) and denote 
the Black-Scholes price of call options at time t by BSt{St, y,T) := StBS{k —log St, y, r) with BS defined in (12.81) . 
Using cr(V) ^ Tt for alH > 0 and (13.11) we see that 


pOO 

dCt{T,k)= dBSt{St,y,T)(p{y)dy 
Jo 


dtBStiSt,y,T) + -S;dssBStiSt,y,T)y ) Cpiy)dydt + St f dsBSt{St,y,T)y/yC,p{y)dydWt 


(B.l) 


= St 


dsBSt[St, y, T)y/yCp{y)dydWt, 


Jo 

where the second line follows from Ito’s Lemma and the last line from the Black-Scholes PDE. Construct now 
a hedging portfolio with a quantity (ft of stocks and ipt of bonds, = 4>tSt + ipt, and therefore (using the 
self-financing condition), since the instantaneous interest rates are null. 


(B.2) 


d-Kt = 4>tdSt = (fty/VStdWt. 


It is not possible to match dynamics between (IB.1|) and (IB.21) since (ft needs to be adapted. Any position in 
the stock will leave a residual exposure for the net portfolio. One could choose a delta position that minimises 
the variance of the hedged portfolio, but we leave this for future research. If one adds a call option C.{T, k) in 
the hedging portfolio tt, with T > t, then 


Tt — (ftSt + 9tCt{T, k) + ift, 
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and C't(T, k) can be perfectly hedged by choosing (jjt = 0, i/’t = Ct/r, k) — 9tCt(T, k) and 

^ _ /o°°5sBSt(5't,y,T)^Cp(;/)di/ 

/“9.BSt(^t,2/,T)VyCp(2/)di/' 

Interestingly one would not hold any position in the stock. 
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